Electromechanical noise in a diffusive conductor 
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Electrons moving in a conductor can transfer momentum to the lattice via collisions with impu- 
rities and boundaries, giving rise to a fluctuating mechanical stress tensor. The root-mean-squared 
momentum transfer per scattering event in a disordered metal (of dimension L greater than the 
mean free path I and screening length ^) is found to be reduced below the Fermi momentum by 
a factor of order l/L for shear fluctuations and {^/L)^ for pressure fluctuations. The excitation of 
an elastic bending mode by the shear fluctuations is estimated to faU within current experimental 
sensitivity for a nanomechanical oscillator. 
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Impressive advances in the fabrication of miniature 
mechanical oscillators provide new opportunities for re- 
search in mesoscopic physics The coupling of elec- 
trical and mechanical degrees of freedom is of particular 
interest. We mention the observation of thermal vibra- 
tion and acoustoelectric effects [Q in carbon nano- 
tubes, the coupling of the centre-of-mass motion of Cgo 
molecules and single-electron hopping and also the- 
oretical work Q on the coupling between a tunneling 
electrical current and a localized phonon mode. 

The present paper was motivated by a question posed 
to us by M. Roukes: Electrons in a metal collide with 
impurities and thereby exert a fluctuating force on the 
lattice. In equilibrium this electromechanical force can 
not be distinguished from other sources of thermal noise. 
Might it be measurable out of equilibrium by driving a 
current through a nanoscale oscillator? To address this 
question one has to consider a delicate balance of forces. 

We will provide both a general theory and a specific 
application to the electromechanical excitation of a bend- 
ing mode in the geometry of Fig. |l|: A thin elastic beam 
connecting two massive Ohmic contacts. The beam could 
be a conductor or an insulator covered with a metal {e.g. 
a metallized suspended silicon beam [Q). We calculate 
the excess noise in the bending mode that arises in the 
presence of a dc voltage V and conclude that it should 
be observable on the background of the thermal noise. 

Let us first discuss the order of magnitude. The noise 
at low temperatures is due to the AfeV/Ep "noisy" elec- 
trons within a range eV of the Fermi energy (with J\f 
the total electron number in the metal). Each electron 
transfers to the lattice a typical momentum Ap ~ in a 
scattering time r. The mean squared momentum transfer 
in a time t for uncorrelated increments Ap would be 



{AfeV/EY){Apf{tlT) ~ Um*eV{tlT) 



.t, (1) 



with m* the electron effective mass. 

We find that Vmax overestimates the fluctuations in 
the transverse force. The actual noise power is of order 



V ~ {I / L)'^'P-max, with L the length of the beam and I the 
mean free path in the metal. The reduction is due to the 
fact that subsequent momentum transfers are strongly 
correlated, basically because an electron being scattered 
back and forth alternatingly transfers positive and nega- 
tive momentum to the lattice. The factor (l/L)^ reduces 
the noise substantially, but we estimate that it should be 
observable in an oscillator with a 10~^^N/-\/Hz sensitiv- 
ity ffl. 
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FIG. 1. Sketch of an elastic beam clamped at both ends 
to a contact and covered by a metal layer. A current flowing 
through the metal excites a bending mode u{x) of the beam. 



We combine two independently developed theoretical 
frameworks: The dynamic theory of elasticity and 
the kinetic theory of fluctuations [|oj. We start from 
the Boltzmann-Langevin equation of Kogan and Shul- 
man. This is a kinetic equation with a fluctuating source 
(5J(r,p,t), that describes fluctuations in time of the dis- 
tribution function n(r,p,t). 



(Sf + V • Vr + eE • Vp + S)n = 6 J. 



(2) 



Here p = m*v is the quasimomentum and E(r, t) is the 
electric field. The collision integral S for elastic scatter- 
ing on impurities (with rate W) is given by 
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5n(p) = (H-(p.p')Kp)-n(p')])p' 



(3) 



The angular brackets indicate an average over the direc- 
tion p' of the momentum p', with |p'| = |p|. 

The noise source 5 J has zero time average and variance 



(5 J(r, p, t)5 J(r', p', t') = (5(r - Y')5{t - t')5{e - e')^-^ 
X [47r(5(p - p')(VF(p • p")(n + n" - 2nn"))^» 
-W{-p-t,'){n + n' -2nn')\. (4) 

Here v{£) is the density of states (at energy e = /2m*) 
and n is the time-averaged distribution. We have set 
Planck's constant /i = 1, so that n is dimensionless, and 
abbreviated n' = n(r,p',t), n" = n{r,p",t). 

The force density f(r, t) exerted by the electrons on 
the lattice is the divergence of a symmetric tensor 11 that 
can be decomposed into an isotropic pressure Hq and a 
traceless shear tensor S: 



(5) 



In the approximation of a deformation-independent ef- 
fective mass one has [ppjl 



a/3 



m* / dpVaVpn. 



(6) 



The time-averaged force (||) vanishes, since it contains 
a derivative of the spatially uniform time-averaged dis- 
tribution n. Physically, this results from a cancellation 
between the electrical current drag on impurities (the so- 
called "wind force" ) and the electric field force exerted on 
the ions Moreover, since f is a total derivative the 
net fluctuating force vanishes as well at low frequencies 
(ignoring boundary contributions). Although the center 
of mass does not move, there are fluctuating compres- 
sion modes (driven by Ho) as well as torsion and bending 
modes (driven by S). The driving force J-{t) for each of 
these modes is obtained by weighing f (r, t) with a sen- 
sitivity function g(r) proportional to the displacement 
field of the mode, 



drig 



dr HoV • g 



dra 



(7) 



The two contributions Hq and S can be separated by 
expanding 7i(r,p,<) in spherical harmonics n'^*)(r,e,t) 
with respect to the direction p of the momentum. It is 
convenient to write the spherical harmonics in cartesian 
(rather than spherical) coordinates. 
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(8) 



Here p*^ is the traceless part of the symmetric tensor 
Pai ■ ■ 'Paq- These tensors form an orthonormal set |]l3t, 



(2m -HI)!! 



(9) 



The tensor A^^™^ projects onto the traceless symmetric 
part of a tensor of rank to. Wc will need A^"'^ — Sap and 

^apa'P' ~ ^^aa'S0f3' + 5(5q/3'<^/3q' " ^SaffSa'/S'- (10) 

In view of the orthogonality of different spherical har- 
monics, one has 

Ho = de2ei^n^"\ ^a/s = ^ J de2evn^^l. (11) 

The two harmonics n^"^ and n*^^) have to be found from 
the kinetic equation (|^). We first consider the harmonic 
n^^^ that determines the shear tensor S, and then discuss 
the harmonic n*-"^ and resulting pressure Hq. 

To obtain an equation for n^^^ wc multiply both sides 

of Eq. (H) by pp and perform an angular average. Em- 
ploying the diffusion approximation on length and time 
scales larger than I and r, we neglect the derivatives with 
respect to t and r in Eq. (||). Also, in the linear response 
approximation, we neglect the derivative with respect to 
p, since it gives a term bilinear in E and j. What remains 
is a local relation between n'^' and the second harmonic 
J(^' of the fluctuating source. 



(12) 



The momentum transport time T2 is defined by 

l/r2 = lj\w{m~e)- (13) 

For anisotropic scattering the time T2 is larger than 
the charge transport time r, defined by 1/r = 
\ J-i d£. W{£_){1 — ^). (For isotropic scattering T2 = t = 
1/W.) 

The correlator of (5J'^' follows in the same way from 
Eq. (^). In the diffusion approximation we replace n, n' , 
n" by n^^\ Using Eq. (0) we arrive at 



<SjS(r,e,t)Wi%,(r',e',i') = AS„,,,5(t - i') 



X dir-r')S{e^e')n'^'^^r,e)[l-n^°\r,e)]. (14) 

Since n^'^)(r, e) differs from or 1 only in a narrow range 
near the Fermi level, we ignore the energy dependence of 
v and T2 and evaluate them at e = Ep. 

We quantify the transverse momentum noise through 
the correlator of the shear tensor, 

/oc 
dtj:ad^,o)j:a'f^'ir',t). (15) 
-00 

Combining Eqs. ([ll|), (12), and (14), we obtain the result 
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C(2)(r,r') = ^{m*vl)^T2iyS{r~r')Kir)A^^\ (16) 
K{r)^ f den<^'^\r,e)[l-n^°\r,e)]. (17) 



The kernel K is given by K = {eV/L)x{l - x/L), 
where V is the voltage applied between the two con- 
tacts at a; = 0, L. The parabolic profile K{x) requires 
k-gT <^ eV and the absence of inelastic scattering. 

We now turn to the pressure fluctuations. Instead of 
Eqs. (^2|) and (|lj) we have the fluctuating drift-diffusion 
equation |14 



j + DVp - crE = er y de z^wW^^^ = SI, (18) 

dIc.{r,t)SIp{r', t') = 2a5^p6{Y - r')S{t - t')K{r), (19) 

which relates the fluctuations in the charge density p — 
e J devn'^^^ and the current density j = J d£vvxi^^\ 
Once we know the charge density fluctuations we can 
find the fluctuating pressure from 



(20) 



cf. Eq. (n). We have introduced the diffusion constant 
D — ^WpT^ the conductivity a = e^vD, and the mobility 
p = er/m*. 

The correlator C'"-* of the pressure fluctuations is de- 
fined as in Eq. (p^, with S replaced by Hq. To close the 
problem we need the continuity equation, 9p/i9i-|- V-j = 
0, and the Poisson equation kV • E = 5/9 (with k, the 
dielectric constant) . The time derivative of p in the con- 
tinuity equation may be omitted in the low-frequency 
regime of interest. The fluctuations in the electron den- 
sity then obey 



D\l'^5p- {cr/n)5p = V -81. 



(21) 



The current fluctuations create a fluctuating charge 
dipole, that is screened over a length ^ = [kD/oY^'^ — 
(nj e^vf-^"^ ■ On length scales ^ ^ one may neglect the 
diffusion term in Eq. ( |2l| ) and use the local relation ||l^] 
8p = -(K/cr)V • 51. Eqs. dH) and @i then yield 



C("Hr,r') 



4(7^4 Q Q 

p? dr dr 



jS{r~r')K{r). 



(22) 



The next step is to use the results ( |16[ ) and 
( p^ ) to estimate the low-frequency noise power V — 
2 /^di J^(0)J^(t) of the fluctuating force T{t) that 
drives a particular oscillator mode. To that end the cor- 
relator (|l^) is integrated over r and r', weighted by the 
sensitivity function of the mode as in Eq. (|^). For a bend- 
ing mode we use Eq. (|l|), which gives the noise power 



p(2)^A(^*,^)V2. /rfrif(r)A^^^„, 



(2) dgp dgp' 



dra dra' 



(23) 



For a compression mode in a metal of size 3> ^ we use 
Eq. (p^) and find 



-pio) 



J drK{r) 



VV-g| 



(24) 



For an order of magnitude estimate, we take K ~ eV, 
g ~ 1, and also estimate spatial derivatives by factors 
1/L and the volume integral by a factor V. For sim- 
plicity wc assume isotropic impurity scattering, so that 
T2 = T. Then the noise power due to fiuctuations in 
the shear tensor is of order T-"*-^-* ~ {m*v^)'^Ti'eVVL~^ 
and the noise power due to pressure fluctuations is of 
order T'^"^ ~ aS^'^ p^^eVVL^"^ . It is instructive to write 
these two estimates in the same form, using the identi- 
ties a/ep — ^m*Vpi' ~ JV /V = fic, with Uc the electron 
density. One flnds 



(25) 



with T^max = Mm* eV /t the noise power for independent 
momentum transfers mentioned in the introduction. 

The experimental observation of the shear tensor fluc- 
tuations looks more promising than of the pressure fluc- 
tuations, firstly because ^ is typically ^ {IL)^/'^ so that 
'pi'i) <^ 7-'(2), and secondly because a typical oscillator 
operates in a bending or torsion mode rather than in a 
compression mode. For that reason we will now limit 
the more quantitative calculation to P^^-*. We consider 
a bending mode u{x) cos ujot in the geometry of Fig. 
The sensitivity function g{x) = u{x)/u{xq) equals the 
displacement (in the j/-direction) normalized by the value 
at a reference point Xq. We choose xq — L/2, so that T 
is equivalent to a point force at the beam's center. Eq. 
(E3) now takes the form 



d-r 

V = ln,pAlA/L) ^Kix)[Lg'{x)r, 



(26) 



with A = V/L the cross-sectional area of the metal layer. 

The wave equation for transverse waves is biharmonic, 
d'^u/dx^ — k^u. The solution for doubly clamped bound- 
ary conditions is 

u{x) — (sin /cL — sinh fcL) (cos fca; — cosh fcx) 

— (cos kL — cosh kL) (sin kx — sinh kx) , (27) 

with the resonance condition cos fcL cosh fcL = 1. We 
use the lowest resonance at kL = 4.73. Substituting 
K = {eV/L)x{l — x/L) and integrating we obtain the 
excess noise V = ^ncPF{lA/L) x 0.83 eF. If we insert 
values typical for a metal, Uc = 10^^ m"'^, pp = lO"^"* Ns, 
I ~ 100 nm, and choose typical dimensions A/ L = 10 nm, 
then the force spectral density at = 1 mV is V = 
10~32 ]\j2^jj2. This is well above the thermal noise power 
at low temperatures (of order 10^"^"' N^/Hz at T = IK 



§)• 



It is instructive to apply the result (|2^) to a sys- 
tem in thermal equilibrium, when K(x) = k^T for all 
X. In this case an independent estimate of the noise 
Vq is provided by the fiuctuation-dissipation theorem: 
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Vo = 4kBTMu!o/Qo, with M the active mass of the os- 
cillator and 1/Qo the electromechanical contribution to 
the inverse quality factor. Eq. (E9) gives 



1 



Nr. 



dx 



(28) 



This electromechanical quality factor might be measur- 
able in a superconducting metal, as an increase in the 
overall quality factor when T drops below the critical 
temperature. One can also calculate Qo directly as an 
"absorption of ultrasound" by conduction electrons 
providing a consistency check on our analysis. However, 
the nonequilibrium noise ( p6| ) with an x-dependent kernel 
K{x) can not be obtained from acoustic dissipation. 

Before concluding we mention an altogether different 
mechanism for electromechanical noise, which is the cou- 
pling of a fluctuating surface charge 5q{t) on the metal to 
the electromagnetic environment. In the presence of an 
electric field Eq between the metal surface and the sub- 
strate {e.g. due to a mismatch in work functions), the 
charge fluctuations will give rise to a fluctuating trans- 
verse force with noise power 

Here C is the capacitance to the ground and R — L/A<t 
the resistance of the metal. The ratio Vcm/V^'^^ ~ 
[EqC L'^ / elMY is quite small for typical parameter val- 
ues. The reason is that the environmental charge fluc- 
tuations are a surface effect, while the whole bulk of the 
metal contributes to T''^^^. Although the noise per elec- 
tron was found to be small in l/L, the total noise power 
■p^^' is big due to the large number MeV/ E-p of contribut- 
ing electrons. 

In summary, we have addressed the fundamental ques- 
tion of the excitation of an elastic mode in a disordered 
metal out of equilibrium, as a result of the fluctuating 
momentum that an electrical current transfers to the lat- 
tice. The effect is small but measurable. The charac- 
teristic linear dependence of the electromechanical noise 
on the applied voltage should distinguish it from other 
sources of noise. We believe that a measurement is not 
only feasible but worth performing. Indeed, the nonequi- 
librium electric current noise has proven to be a remark- 
ably powerful tool in the study of transport properties 
I prt , precisely because it contains information that is not 
constrained by the fluctuation-dissipation theorem. The 
noise considered here could play a similar role for me- 
chanical properties. 
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